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Abstract 

In this paper we introduce a diagrammatic equation for the planar 
sector of square non hermitian random matrix models strongly remi- 
niscent of Polchinski's equation in quantum field theory. Our funda- 
mental equation is first obtained by a graph counting argument and 
subsequently derived independently by a precise saddle point analysis 
of the corresponding random matrix integral. We solve the equation 
perturbatively for a generic model and conclude by exhibiting two 
duality properties of the perturbative solution. 

1 Introduction 

Random matrix models have wide applications in subjects ranging from sta- 
tistical physics PP to two dimensional quantum gravity [2]. They have long 
been related to triangulations of two dimensional surfaces pj and various 
combinatorial counting problems [ll|5l|6j have been addressed with their help. 
Techniques inspired by quantum field theory have been adapted to such mod- 
els [7] , and non identically distributed random matrix models proved crucial 
in the study of non commutative quantum field theories [51 El [IHl E] • 

For a large class of random matrix models, in the limit of large matrices 
(the so called "large limit"), the behavior of the partition function (and 
all correlation functions) is dominated by a saddle point. The planar sector 
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[3] is the first contribution to the evaluation of the random matrix integral 
by the saddle point method. The study of this saddle point is best addressed 
by means of the resolvent function which obeys a certain quadratic equation. 
In the large limit the resolvent is the generating function of planar graphs 
with one vertex, thus having a very natural graphical interpretation. 

Rectangular non hermitian random matrices been extensively studied (see 
[5] and references therein). However the treatment presented [5] is not imme- 
diately applicable to square matrices: it is singular already for a free model 
(as it will be explained in detail in section |4]). The solution implemented 
in [3] is to study rectangular non hermitian matrices and take the limit of 
square matrices after having taken the large N limit. 

In this paper we propose a new approach to square non hermitian ran- 
dom matrices. Our approach does not suffer from the singularities of |5] 
(and sheds light on there origin). A detailed analysis allows us to derive 
a tower of equations obeyed by the coefficients of resolvent of such models 
both by a graph counting argument (inspired by the Polchinski equation) 
and by a saddle point analysis. We will solve these equations in perturba- 
tions and prove certain duality properties of the solution. Our results encode 
the census of planar diagrams with one vertex generated by an arbitrary non 
hermitian random matrix model (which can be translated into the census of 
planar connected amputated graphs with one external face). Non hermitian 
matrix models have been studied in [T^ |T2] by the method of the "loop in- 
sertion equation" (very similar to our equation for the resolvent) which can 
be derived either by a combinatorial argument (for the resolvent itself) or 
by a functional manipulation. In [T2l [T3] the equation is solved by contour 
integrations (allowing rapid access to large order behavior of the solution) 
rather than trough the "brute force" approach we take here. Their results 
provide a check of our method and suggest possible roads to generalize our 
tower of equations for non planar graphs. 

In higher dimmensions, random matrix models generalize to group field 
theories [El |T3], associated to cellular complexes [IB] with boundary [T7] . 
A generalized notion of planarity adapted to group field theory has been 
proposed p!8], and various power counting results [HI |19l EQ] have been 
established. Our long term goal would be to generalize the results of this 
paper for such theories. 

This paper is organized as follows. In section [2] we present the non hermi- 
tian random matrix models and the relation between the resolvent and the 
connected planar functions with one external face. In section |3] we derive 
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our fundamental equations by a graph counting argument, and recover them 
in section |4] by a saddle point analysis. Section [5] provides the perturbative 
solution of these equations. In section [6] we exhibit two duality properties of 
the perturbative solution. Section [7] draws the conclusions of our work. 



2 Non Hermitian Random Matrix Models 

In this section we review in some detail the non hermitian random matrix 
models. We use this opportunity to introduce our notations, starting with 
the most important: troughout this paper we will denote 

^ p\{n-p)\ ' 

the binomial coefficient. 

Let M be an square (A^ x N) non hermitian (M ^ M^) matrix. A non 
hermitian random matrix model is a probability measure 



df^ViaM^M) = [^J [11 ^^^b dMab 

ab 



^Tr(MtM)+Tr[V(oMtM)]^ 



oo 

V{aM^M) = -J2- 9p(^^ {M^MY . (2) 



1 P 



with a some small coupling constanlj^ The normalization is chosen such that 
at a = 0, d^iQ is a normalized gaussian measure. The partition function and 
the correlations of this measure are denoted 

^ = (1) = J df^V(aMlM) 

Ml . . . M,, . (3) 

The arguments of a correlation, . . . M^^, are called "external points" . 
The partition function and the correlations are evaluated as sums of oriented 
ribbon graphs with external points build as follows 



^The exponent Tr(A/'fM) + Tr[F(aM^M)] in eq. ^ is called the action of the model. 
A rescaling of both M and by brings it in the more familiar form iTr(AfUf) + 
Ti[V{M^M)]. 
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Figure 1: A typical ribbon graph. 



• Vertices (drawn as crossroads of ribbons) are generated by VlaM'^M). 
Each vertex has 2p alternating halflines (out of which p are Af'^'s and 
p are M 's) and a weight A^^ gp a^. 

• Lines (drawn as ribbons) fall in two categories: internal and external. 
The internal lines connect a halfline and a M halfline. The exter- 
nal lines connect an internal M"^ (or M) halfline with an external M 
(respectively M'^) pointj^ A lines has a natural orientation (say from 
Mt to M) and a weight ^. 

• Faces are closed circuits formed by the sides of the ribbons. They also 
fall in two categories, internal and external. The internal faces are 
closed circuits which do not pass trough any external point. They have 
weight A^. The external faces (containing external points) have weight 
1. 

The identification of the external faces of a graph is slightly non trivial. 
To count the external faces of a graph one "pinches" the external points, that 
is one connects the two sides of the ribbon arriving at any external point. 

Figure [TJrepresents a typical ribbon graph. The "pinching" of the external 
points is represented by the dotted lines. This graph has two external points, 
two vertices, three internal lines, two external lines two internal faces and 
one external face. 

The partition function Z = (1) is the sum of graphs with no external 
points (vacuum graphs). The non trivial correlations have the same number 
of M and M"^ external points. They write 

^AmUK',.'. ■ ■ ■ Ml^.^M,^^,^) = , (4) 

\ ' g 

^We consider the lines connecting two external points as external. 
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where Q are all ribbon graphs with no vacuum connected components (due 
to the division by (1)) and A{Q) is the amplitude of Q. 

To compute A{Q) we denote V the number of vertices (of coordinations 
2pi, . . . 2py), L the number of internal lines, the number of external points, 
F the total number of faces and B the number of external faces of Q. We 
have the two topological relations 

2^p^- Ne = 2L V-L + F = 2-2g, (5) 

with g the genus of the graph Q (thus for the graph of figure ^ g = 0). 
Furthermore, for an external face b of Q, we index the Nf, external points M"^ 
by i^''^ and the external points M by We associate to the face b the 
ordered set 

and we have 

6i n 62 = V61 62 

...^%} = {l...N,} [jiA"^ ...jSI} = {1... iVj . (7) 

b b 

With these notations the amplitude of Q writes 

A{g) = K{g)a^^P^N^-^'^+''^^+^^-^^f[gp^ 

v=l 

where K{Q) is some combinatorial coefficient. 

Of particular interest in the sequel is the resolvent defined as 

""W = ]v (T) °j^„^'N « ■ 



We will denote its limit when — > 00 by uj{z). The coefficient 

1 1 

Nil) 



dpm = ^r—A^M^ Mf]) , (10) 
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is the sum over all connected vacuum graphs having a special vertex of co- 
ordination 2p and weight 1 (the trace and matrix product over the external 
points allows us to reinterpret them as the halflines of this special vertex). 
The global power in of a graph contributing to (ip;(Ar) is, using eq. ([8|, 

i_^(y-l)-L+F ^ ^-2, ^ 

and we recognize the classical result that dp = limN^oo dp- {n) is the sum 
over connected planar vacuum graphs with a special vertex of weight 1 and 
coordination 2p. 

Before concluding this section we present the relation between the con- 
nected planar one vertex functions dp and the connected amputated planar 
functions with one external face. Our motivation is twofold: first, from a 
quatum field theory perspective, the latter dominate the floating Wilsonian 
actior][^21j and second, the reasoning we present below is similar to the one 
we will use in section |3] to derive our fundamental equation. 

A planar connected amputated function with one external face and 2t 
external points writes 

N''(mI,M,,^, . . . M,,^,y^''''^' = Gt 6,^y^6^,^^ . . . . (12) 

Consider a graph Q contributing to dp (depicted in figure |2]). ^ is a 
connected planar vacuum graph with a special vertex of coordination 2p 
and weight 1. The special vertex is represented by dashed lines in figure |2] 
Consider a halfiine of the special vertex, denoted Ai in figure [2j As ^ is 
planar, if we erase the special vertex Ai becomes an external point of some 
planar connected component C with one external face, represented as shaded 
in figure |2} Say the C has 2t external points. All of them originate from 
halflines on the special vertex, Ai, A2, . . . A2t- We denote 2ki the number of 
halflines on the special vertex between Ai and A2, 2^2 the number of halflines 
between A2 and ^3 etc. up to 2k2t the number of halflines between A2t and 
Ai. As Q is planar, the halflines 2ki must connect among themselves and 
therefore form a graph corresponding to dk^. The same is true for 2/c2, 2^3 

■^To be distingushed from the effective action which is generated by one particle irre- 
ducible graphs. 
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etc. Therefore dp writes 



k2t 



(13) 



t=\ 



ki+...k2t=p-t 




Figure 2: A graph Q contributing to dp. 



Equations (13) can be solved iteratively for Gt- The first equations are 



di = Gi 

(^2 ~ G2 ~\~ '2Gidi 

ds = G3 + 4:G2di + Gi{dl + 2d2 



(14) 



which are solved by 

di 



Gi 

G2 ~ ^2 — ^Gidi = — 2(i^ 

G-i = d-i-4:G2di-Gi{dl + 2d2) = d3-6did2 + 7dl 



(15) 



Before concluding this section we introduce a notation. As V{x) = 
— p9p^^ is a formal power series, we will denote its derivative (in the 
sense of formal power series) by 



(16) 
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3 The Diagrammatic Equation 



In quantum field theory the Polchinski equation [22j is an equation for the 
floating Wilsonian action S at scale A and writes 



OaS 



6S 



-dAK{x,y) 



6S 

mi) 



6^S 



6(f){x)6(f){y) 



dAK{x,y), (17) 



where K{x, y) is the propagator with ultraviolet cutoff A. The action S is 
the sum of connected amputated graphs, and the propagator K is associated 
to lines. This equation naturally translates in a graph classification: a line 
in a graph contributing to the action S can either be a "tree line" separating 



two distinct effective vertices (the first term in eq. (17)) or it can be a "loop 
line" for some effective vertex (the second term in eq. (17)). Adapting this 



idea to ribbon graphs will lead to our fundamental equation. There are two 
important aspects one needs to take into account. First, in stead of dealing 
with the floating action, we will deal with the resolvent of the matrix model. 
Second, we will look for an equation valid in the N ^ oo limit, thus involving 
only planar graphs. 




special vertex 



2q-2k halfiines 

Figure 3: A loop line on the special vertex. 

First consider a free model, obtained by setting a = in eq. We will 
denote (the — oo limit of) the resolvent of the free model co^'^'lz) and its 
coefficients d^pK 

With this notations, d^^Ji is the sum of all planar vacuum graphs with one 
special vertex of coordination 2g + 2 and weight 1. Consider a line starting 
on the special vertex. As the graph contains no other vertices, this line must 
also end on the special vertex. It is therefore allways a loop line for the 
special vertex. The line encloses a certain number of halfiines, say 2k, and 
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leaves on the exterior 2q — 2k halflines as in figure |3} As the initial graph is 
planar, the interior (exterior) halflines necessarily connect into planar graphs 
with one special vertex of coordination 2k (respectively 2q — 2k). Therefore 



d 



(0) 



fc=0 



(0) ,{0) 
k "-n-k^ 



d'. 



(0) 



(18) 



The attentive reader will note that up to now we did nothing new: eq. (18) 
is nothing but equation eq. (13) supplemented by the condition that the 
connected amputated functions of a free model are trivial, Gt = Sn- 




special vertex 



Figure 4: A tree line in the interacting model. 

For an interacting model, a line originating on the special vertex of dg^i 
can either (as for the free model) end on the special vertex, or it can end 
on a gn vertex as in flgure |4| Up to a weight factor, the graphs contributing 
to the latter are in one to one correspondence with the connected planar 
vacuum graphs with a special vertex of coordination 2q + 2n. To compute 
the mismatch in weight factors, recall that the weight of the vertex Qn is 
^gn.Ci"' and note that the line has a choice of one among the n halflines of 
appropriate orientation on Qn, therefore 

q oo 

dq+1 = ^ dkdq-k + ^ QnOl^dq+n ■ (19) 
k=0 n=l 

This are our fundamental diagrammatic equations. The rest of this paper 
is devoted to their analysis. 

4 The Saddle Point Analysis 

In this section we derive the equations ( [l9| by an independent method. First, 
by a carefull saddle point analysis adapted to a square non hermitian random 
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matrix model, we will derive an equation obeyed by the resolvent u}{z) and 



then, by developping in powers of z, we will rederive the equations (19) for 



the coefficients dp. 

To evaluate an integral with measure ([2]) by a saddle point method one 
needs to change variables from M*^, M^^ to the eigenvalues of M'l'M |23]- 
At this stage one has a choice. One can either chose to change variables to 
Xi the eigenvalues of M^'M, or to change variables to the eigenvalues of 
y/ M'^M. For example in [5] the author chooses the former. However, for 
square matrices, there is a heavy price to be paid in doing this choice. 

To understand the problem consider first Euler's gamma function 

POO 

T{n) = / r-^e-'dt . (20) 
Jo 

It is well known that this integral is governed by a saddle point to = n — 1, 
and the saddle point evaluation of it proves Stirling's formula. However a 
polynomial change of variables u = yields 

T{n) = - / e-""rfM , (21) 
^ Jo 

and in this new form the integral is not governed by a saddle point: the 
Jacobian of the polynomial change of variables absorbed it. 

This same phenomenon occurs for square non hermitian matrices. Fol- 
lowing [3] one writes (up to some normalization) the partition function 



•^0 7. • - • 



(22) 



k i<j 



For simplicity let us set a = 0. Raising the Van Der Monde determinant in 
the exponent yields 

(23) 

i j<i 



The intergal (22) is not governed by a saddle point: the would be saddle 
point equations 
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are incompatible: adding them all up yields = 0. Considering, as in 
[5], rectangular matrices solves this problem. In that case the saddle point 
equations pick up an extra term and are compatible. As previously stated, 
their results can be applied for square matrices only after taking the large A^ 
limit. 

In the sequel we will take the second alternative and use the variables 
Hi, eigenvalues of V M^M. The Jacobian of this change of variables has 
been computed in [23] and the partition function writes (again up to some 
normalization) 



k i<j 

-NF{{y,}) 



Yldyf 



k e 



where we denoted 

Fi{y^}) = T.y" + Il ^("^') - ^ E 1- - 1^ E 1- \y" - y' 



N 

We have the following lemma 
Lemma 1. uj{z) obeys the quadratic equation 

u{zy - P{z)uj{z) + Q{z) = 

with 

P{z) = 1 + aV'{az) 
Q{z 



(25) 



(26) 



i<j 



(27) 



1 a 



zV'{az) — mfV'{amf) 



z — mj 



(28) 



Proof: The proof is a variation on the standard derivation of the quadratic 
equation obeyed by the resolvent of a hermitian matrix model. 
We start by writing the resolvent as 

1 1 /„ / 1 



ujn{z) 



N{1) 



Tr 



z-um 



1 


Jo 






1 (,-NF(,y,) 

z-yf 


N 


ir Uk dvk 





(29) 
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By slight abuse of notations we denote the leading contribution (in powers 
of j^) in eq. (29) also by ujNiz). At large both the numerator and the 
denominator above are dominated by the same saddle point yi = rrii, with 
rrii the solution of the equations 



oVi N yi 



2 2yi 



N 







(30) 



where V is the formal power series in eq. (16). The resolvent (at leading 
order in j^) is then 



un(z) = — — 

^ ^ N ^ z-m 



2 ■ 



(31) 



Eq. (25) ensures that all rrii are strictly positive. In order to find the 
equation for ujn{z) we note that any linear combination of the saddle point 
equations is zero, hence 



d 



^ _ 1 ^ 1 1 

^ z — mf 2mi V dy, 

T7 — fl + aV^arri^i 
N ^ z ~mj\ 



yi=mi 



(32) 



1 1 2 

2N^ ~ N 



^ mf — m?- 



In the case of a hermitian matrix model one employs a similar trick, but 
without dividing by 2mj. Each term in the last line of eq. (32) receives a 
different treatment. The first term (also present for hermitian matrices) is 
identified as 



A^ z — mt 



(33) 



The second term (again present for hermitian matrices) rewrites by adding 
and subtracting aV'{az) as 



— V — ^ (aV'iamf) - aV'(az) + aV'(az)) 

i ' 

\ /X 1 ay (am^) — a\^'(a;z) 
= aV\az)uj^{z)^-Y,— 



z — m; 



(34) 
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The third term receives a similar treatment namely it is rewritten as 
1^1/1 1 1 



z - y- \y 



z z 



^ i-M-) + ^ i-«(0) . (35) 



This is a new term, not present in the hermitian matrices case. Finally the 
last term computes (as for of hermitian matrices) to 



Ar2 



EE 



11 1 
mf — m - z — mf 



1 dUN^Z 

N dz 



(36) 



Substituting everything into eq. (|32|) yields 




UNiz) + aV'iaz)uNiz) + ^2^— ' ^ 



z — 

1 duNiz) 



(37) 



The new terms of equation (35) prevent us from taking — )■ oo at this 
stage: although Co'Ar(O) is well defined for any finite A^, its limit ci;(0) is ill 
defined hence one can not conclude anything about the limit of ^a;Ar(0). To 
get around this problem we use again the saddle point equations ( 30 ) to write 







(38) 



Substituting eq. (|38|) into eq. (|37|) yields 
= 



1 duN(z 



1 1 



dz 

a^V'{ayf)-V'{az) 
AT 2-^ 



^N{zf + un{z) - i^~^n{z) + aV\az)u}N{z) 



z — 



a 
N 



Ef'i 



am. 



1 



(39) 



which rewrites after rearranging the terms as 

1 dujN{z) 







dz 

a 

'n 



1 1 

2iV^ 



a ^ f V gmf - V (az) 1 ^ . 

N ^\ z-nij z ' 



(40) 
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In the large N limit the first two terms vanish and a straightforeward com- 
putation proves the lemma [T] □ 



Recall that u{z) = -^dp, therefore lemma [T] can be translated into 
an equation for the coefficients dp. In the reminder of this section we will 
show that this equation is exactly our diagrammatic equation. 



The quadratic equation (27) for ijj{z) can be written as 

1 a sr^ zV'(az) — mJV'iamf) 



uj{z) ^'^{z) = —aV' {az)u{z 



zN 



E 



z — m; 



(41) 



and substituting uj{z) = Xlp=o^ ^ ^^p^ ^^^^ hand side becomes 

oo ^ ^ oo oo ^ 

p=0 p=0 (?=0 



(42) 



The right hands side in eq. (41) needs a bit more work. We first express 



the coefficients dp in terms of the saddle point solution mi as 



°° 1 1 1 

p=0 



dn 



z — mj 



(43) 



and substituting V from eq. (16), the second term on the right hand side of 
(41) becomes 

1 a 



zN 



i n=l 

oo n—l 



z — mj 



^ oo n—l 



n=l 



t=0 



,n-l-t}_ 

N 



2\t 



(44) 



n=l 



t=0 
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Using again the definition of V, the riglit liand side of (41) becomes 



n-l 



a 



t 

^+1 



n=l 

oo 



t=0 
oo 



n=l 

oo 



i=0 

oo 



n=l t=n n=l q=0 



+1 



oo oo 
q=0 n=l 



(45) 



Equating the coefficients of 2; ^ in eq. (42) and (45) yields 

q 00 



(46) 



k=0 



n=l 



which is exactly our diagrammatic equations ( 19 ) 



A selfconsistency check is provided by the the free model (a = 0). In this 



case eq. (27) writes 



(47) 



as uiz 



^00 z ^. The eigenvalues distribution of the free model reads 
P(A^) = ^HX' + ^0) - co{X' - .0)] = ^J2 ^(^' - 



2m' 



and we recover Wigner's semicircle law 

p{X)d\ = p{\^)d{\^) = -V4^A2 dX . 



(48) 



(49) 



5 The Perturbative Solution 

In this section we will solve our diagrammatic equation in perturbations. We 
start by the free model, and subsequently work our way to the solution of a 
generic model. 
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Lemma 2. The solution d^^ of the free model (equations (IS)) is 



(0) _ _J_fy2k W7. 



(50) 



Proof: A direct computation shows that (ig^^ = 1. We proceed by induc- 
tion. Suppose that 



d! 



(0) 



A; + 1 " 



(51) 



Eq. (18) then reads 



" 1 1 

k=0 ^ 

f2^U + l^g-A; + l 

k=0 ^ 

, 2^fc + l ^ 



1 



-2fe 



(52) 



where in the second hne a change of the dummy variable identifies the two 
terms. Separating the term g = in the sum gives 



q + -Z\q + 
which further computes to 



1 f 2 2q , \^ 2 ^2fc^2g- 



2-2fc 



(53) 



,(0) _ 1 / 2 2g I 2 2fc^2g~ 



q-1 



2-2k 

-k 



k=0 
2k 



A \^ ^ f^2k \ ^2q-2-2k 

""l^k + l " (q-l)-k + l '^-'-^ 



k=0 
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Taking into account the first and third hnes in eq. (52) allows us to write 

2 



d 



(0) 
<?+! 



q + 2\q + l ^ 
q + 2\q + l « 



+ 4(g + l)4°)-4rff 



q + 1 



(54) 

□ 



We now turn our attention to the interacting model. We will solve the 
equation (19) in perturbations around the a = solution we just derived. 

1 



We substitute 



^9 = 4 



in eq. (19) to get 



q oo 



= 0,i>l 



(55) 



E 4ii"^ = E E ( E 4^^4-7^) «^ + E E ^"4r:^ 



s=0 k=0 s=0 i=0 s=l n=l 

Equating the powers of s on the left hand side and right hand side yields 

it = E E ^iP^ + E 9n47:' ^ 1 • (57) 



k=0 i=0 



n=l 



This equation mixes all orders of perturbation from (0) to (s). Solving it 
comes to write the terms at order (s) as functions of terms at lower orders. 
To this end we define 



q s — 1 



T'q% = E 9ni7:' + E E 4^^4- ^ > 2 , m 

n=l k=0 i=l 

such that at the order (s), T^^'^ depends only on lower orders d^'^^ to 



Lemma 3. The solution of (51) at the order of perturbations (s) is 



9+1 ~ Z-^ ^q-p -^p+i 



(59) 



p=0 
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Proof: The proof of this lemma is somewhat technical. We start by inserting 
the definition of T^^^ in eq. (57) to get (for s > 1) 

45i = ^i:U2x:41*^- (60) 

We will solve this equation as a recurion over q. However note that the 
right hand side involves all d^j^^ with k < q + 1. We will first derive d^^'^i as 

a function of only d[^^ and T^^^ Then we will solve this recursion and find 
dg^li as a function of only T*^'*^ 



Step 1: Recalling that d^ 

1 



for s > 1, we have 

9-1 



fc=l p=0 



Reinserting eq. (60) yields 



9-1 



- T^'^ + 2 V ( T^"^ + 2 V d^^^ d^ 

p=0 t=0 
p=0 p=0 t=0 



t=0 p=t 



fc=i 

9 



fc=l 



t=0 u=0 



The sum over m reproduces the right hand side of eq. (18), thus 

9 9-1 



(0) As) 



k=l 



t=0 



Using dl^^ = 1 and again eq. (|60|), we have 



4t. = r<:>+2$:4vr+4(EW-<') 

/c=l t=0 
9 



As) _ rp{s) 

-fg+l 



4^^ 



(61) 



(62) 



(63) 



(64) 



k=l 
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and regrouping S^li on the left hand side yields 

4i = 44^^-2E41^r+^a- 



(65) 



k=l 



(s) 

Our first objective is now achieved, as dg^i is expressed solely in terms of 
4'^ and T(^). 

Step 2: Reinserting eq. (65) into itself p times we find 



As) 
^q+1 



q-p p 



= 4'«4-',-E(E2'*-<i-.)^-' 

k=l s=0 
q q-k 



''q-k- 



rp{s) I rp{s) 

^k 



(66) 



k=q—p+l 



s=0 



is) is) 

Using d\ =T{ \, setting p = q — 1 and grouping the first two terms yields 



q—k 



k=l s=0 



(s) I rp{s) 

q+1 ■ 



(67) 



By now we have expressed only in terms of T^'^K To simplify this 
expression, we use the result of appendix [A] 



C _ V ^^(0) - If] ^_f^2q+2\ 

^i~2^4k+i^k - 2\ 49+1 ' 

k=0 



(68) 



and write 



4i = E(4^^'"'-2-4^"'^'v.)^f 

k=l 

= Eh'^'"''~2 

fc=l 

and finally 



(69) 



AQ-k+l^f] ^2g-2fc+2 1 



(s) 



g+1 ' 



^(^) _ ^2q-2k+2rp(s) rp(s) _ ^2q-2prp 



is) 

p+1 • 



(70) 



p=0 
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□ 



For practical computations it is useful! to recast the perturbative solution 
of lemma |3] in terms of generating functions. 

5.1 Generating Functions 

Let us define the generating functions at order (s) 



-^"(-^ = E ;^4^^ ' = E • (71) 



q=0 5=0 



where we set Tq = for all s > 1. The solution of the non interacting 
model of lemma |2] gives us the generating function at order zero 



5:^^^.-4(1-^1-1). (72) 

g=0 ^ 

For all s > 1 = 0, thus lemma [s] translates for the generating functions 
at order of perturbations (s) as 



9=0 q=0 q=0 p=0 



p=0 q=p 



The last ingredient one needs is to translate the definition of T^^+i (eq. ( 58 ) ) 
for generating functions. We first treat the case s = 1 



q=0 q=0 

OO _ ^ 

= (74) 



^1 

g=0 
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Then, for all s > 2 we have 

5=0 9=0 n=l k=0 i=l 

n=l 9=0 j=l 9=0 k=0 

s n~l ^ s—l 



g s—l 



n=l 



q=0 



i=l 



The equations (72), (73), (74) and (75) allow one to compute the gen- 
erating functions w^order by order in (s). In appendix [b] we use them to 
compute uj^^\ w*-^^ and uj^^\ 



6 Dualities 

In this last section we will exhibit two duality properties of the perturbative 
solution. In fact dp is a power series not only in a but also in gk, that is it 
writes 



n-fc 



ni=0n2=0 nfc=0 

j(ni+2ri2H \-knt+--) fiyc\ 

and the coefficient d ■■ is the number of unlabeled, planar, ori- 

entable graphs having a special vertex of coordination 2p and rii vertices gi, 
n2 vertices g2 and so on and so fort. 

It is clear intuitively that these coefficients are not independent. First, 
the number of unlabeled graphs with a choice of vertices is a combinatorial 
quantity and does not depend on which of the vertices is the "special vertex" . 
Second, the vertices of weight gi have coordination two thus they can be 
viewed as decorations on the lines of a simpler graph. We present below the 
precise duahty properties of the coefficients d^^^^-^""^nj!''~^"'' mirroring these 

two observations. To this end we start from the definition of dp 

rfp=^^(Tr[(MtM)^]), (77) 
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and denoting (. . .)o the gaussian {a — 0) correlation, dp writes 



where the superscript indicates that only planar graphs contribute to the 
Gaussian average. In the sequel, we only count the planar contributions in 
all Gaussian averages, but, to simplify notations, we drop the superscript. A 
perturbative development in all the coupling constants gk yields 

1 1 °° 

rl — A7-ni+...nfc+... niH \-knk+- 

- N(l) ^ 

ni,...nk,—=0 

ni\\l) "'nk\\k) 

'tt[{M^ Mf]Tr[{M^ M)y'' . . . Tr[(M"fM) . . . , (79) 
which allows us to identify 

„ , s 1 1 ]\jni+...nk+--. 1 

^(ni+2n2+...fcnfc+...) _ i i iV i 



92 ■••3fe ••• iV (1) ni! . . .Ufc! . . . I"i2"2 . . . A;'^* . . . 
X ^Tr[(M"fM)^]Tr[(MtM)]"i . . . Tr[{M^ M)'']'"' . . . . (80) 

In the gaussian average above the special vertex does not play any distin- 
guished role. 



Tr[(MtM)^']Tr[(M^M)]"i . . . Tr[(M^M)*=]"'= . . 

Tr[(MtM)]"^..Tr[(MtM)f]"^+^..Tr[(MtM)'=]"'=... ) . (81) 



The only thing which distinguishes between the special vertex and any other 
vertex is the combinatorial weight. Balancing this combinatorial proves the 
first duality property 

^ ^ r(ni+...pnp+...q{nq+l)+...) _ ^ 1 Ani + ...p{np+l)+...qnq+...) /qqN 

pUp + l p;9i --gp -Qq ■■■ qUq + l q;9i -ap ■■■9q ■■■ 

for all q with nq> 0. 
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To prove the second duality we introduce the modified Gaussian measure 

-xAf (Tr(MtM) 



ab 



(83) 



and we denote the correlations of this measure {■ ■ ■)o;x- With the help of the 
modified measure we can write 

1 1 ]\fni+l+...nk+... 



d 



{ni+l+...knk+... ) 



p;9i 



ni +1 n 



9k-- N {!) {ni + l)\ . . .Ukl . . . I"i+i2"2 . . . A;"* 
Tr[(MtM)P]Tr[(M"fM)]"i+^ . . . Ti[{M^ M)'']'^^ . . . 

1 1 ]\fni+l+...nk+... ^ 



N (1) {ni + 1)\ . . .Ukl . . . I"i+i2"2 . . . fc"* 



1 d 

Ndx\ 



(Tr[(MtM)P]Tr[(M"fM)]"i . . . Tt[{M^ M f]""^ . . . 



0,x 



x=l 



as the derivative acting on the gaussian measure produces an extra Tr(M"''M) 
insertion. But the correlations {■ ■ ■)o;x of the modified gaussian measure are 
sums over the same graphs as the usual gaussian correlations (. . .)o, with 
the only difference that the lines of the graphs have a weight ^ in stead of 
jj. All graphs contributing to the averages in (84) have the same number of 
lines, namely p + Hence 



,(n.+l+....n,+...) ^ 1 1 iV"^ + l+...n.+... 

p;9-^+\..g:''... N{1) (ni + l)!...nfc! 



1 



1 d 

Ndx 



x=l 



Tracking again the weight factors one has 



d 



(ni+l+2n2+--+fcnfc+...) 

Til +1 nt. 

p;9i ■■•9fc ■•• 



(X m no ^i- 



(85) 



(86) 



These dualities are very usefuU in concrete computations. Using the re- 
sults of appendix [Bl dp writes up to third order in a as 



dn 



q + l 



q + l 



+ aV^'-^^Cl'' + a'9,92 3qCl^ + a%^C'^''. (87) 
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The first duality (1 



82 ) implies for instance 



(3) 
2,93 



7(2) 
^3,92 



(3) 



.93 



(1) 
3,91 



2 1,9192 



2^1,92 



■^2,91 ' 



whereas the second duality (86) implies 



9,9192 
9,9? 



i + ^^(2) 

3 



<?,9? 



g + 1 



"g,9i "g,9i 



(89) 



7 Conclusions 

In this paper we introduced a tower of diagrammatic equations for the pla- 
nar sector of a square non hermitian random matrix model. We solved the 
equations in perturbations and proved two duality properties of the solution. 

Our results can be generalized automatically to rectangular non hermitian 
random matrix models. The generalizations to hermitian matrix models is 
a bit more subtle as the latter generate also vertices of odd coordination. 
We are however confindent that one can find corresponding diagrammatic 
equations related to the quadratic equation obeyed by the resolvent. 

Our equation encodes the census of oriented planar diagrams in a very 
transparent way. Although much more involved, a generalization to graphs 
of arbitrary genus should be possible. 

Another direction of future research is to generalize these equations to 
group field theories. In order to do this several major points still require 
clarification. First, although we have a proposition for a generalization of the 
notion of planarity (type 1 graphs in [18]), we still need a proof that only these 
graphs dominate the partition function. Second, in higher dimmensions, one 
deals with random tensors, not with random matrices. In order to perform 
an indepth saddle point approximation one would ideally need to find an 
appropriate generalization of the notion of eigenvalues for tensors. 
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A A simple sum 



Lemma 4. Let Sq = XlLo 4*^4°^- 



en 



1 



S.= 2[^-^.CX.') (90) 



Proof: Recall that d'f^ = ^C^?. We have the Taylor expansion 



1 °° 1 

l(i-vr^) = $:J^4«>.««. (91) 



2 

and therefore the generating function of Sq is 

oo oo 9 



g=0 q=0 k=0 k=0 q=k 



fc=0 k=0 
1 1 



(0)^fc+l 



2\i-x v^r^/ ■ 



(92) 



But 



-L= = l + y ArCo+f ^"^^ , (93) 

* (}=0 



hence (90) holds. □ 
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B Perturbative computations 
B.l First order 



At first order, eq. ([72]), ([73]) and (|74|) yield 

1 1 1 



and using 



1 °° 1 

i ~ 7 9=0 



we have 



uji^){z) = ^ify lief ---V^-c^'? 

g=0 g=0 

oo 



5=0 



Therefore 



= /l (2g + 2)(2g + l) _ X ^ ^ 



B.2 Second order 

At second order we have 



7(1). 



^(0) 



T^'\z) = g^[j'\z)-^)+g,z[j'\z) 

which simplifies to 

T^'\z) = g,J'\z) + g2zU'\z) --- \) + J'\z)J'^ 

V z z'^ J 
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We group together the terms containing gi 
which compute to 



(100) 



9l 



2J1-I zjl ^ 2 



2J1 



z. 1 



1 1 

=^ + 2 



9l 



1 /I 1\2 1 

~ 4J 



thus 



CO 



1-^V2 z 

z 



9l- 



9l- 



(101) 



/111 1 \ 

+ ^=~0 ^= 



(102) 



Using 



(1 - ^ j 



2:9 ^2 - z'- 
t=i 



(103) 



a; 



•-^^(2;) writes 



9l 



00 ^ 



?=1 



^9+1 



9=0 g=0 q=Q 



(104) 



A direct computation shows that the terms z^z'^^z ^ cancel. Shifting each 
sum by the appropriate amount yields 



00 ^ 



9=0 



^9+1 



+sAT.\c^.t'^.-Y.Cf.f^.-T.Cf^) . (105) 

9=0 9=0 9=0 



9=0 
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and we conclude 



9i 



r^2q 

2^^ 



9i 



2''' + 



B.3 Third order 

To the third order we have 

T(^) {z) = g,J'^ (z) + g,z {j'\z) - + g,z' (c.(°) {z) - \ 
We start by computing 



9\ 



9\ 



1/2 



therefore, factoring co^'^^ in the first and last term of (107) yields 

1 



T^'\z) 



9i92 



z 



2^1- 
z^ z' 



z 

^~2 



1 



1 - 



2 - 2Vi-;---i-' 



1 



^91 



+92Z 



1/2 



1 



9l 



1 



2. 1 



1 

^ ~ 2 



zm ^ 

1 



3/2 



z\ 1- 



The coefficient of gig2 is 

z z 



z 

X 



2J\ 



z 

^ ~ 2 



1 

4 Z 
z 

1 



1/2 
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which, after a tedious but straightforeward computation, simphfies to 

6 1 



We can now write co^'^^ as 

,{3) 



9l- 



z 



1 / 2 



(111) 



+ 9192 



+ 93 



6 



11/^41 1 



1 _ 4 v2 2 \/ " - - -2 



(112) 



Using 



1 



5/2 



1-^ 

z 

the coefficient of in (112) is 



- (t + 2)(t+2) 2*+4 1 

2^ - 



t=0 



12 



(113) 



- (t + 2)(t+2) 2*+4 

2^ ^> 



t=0 



12 



1 f^ (t + 2)(t + l) 1 



t=0 



g=2 

9(9 + 2) 2a 1 



q=0 



while the coefficient of (71(72 in ( |112[ ) computes 



^2 4^'/' 



(115) 



q=0 



Finally the coefficient of (73 in (112) 
z'' ^ _ ^ 1 



2 1 



2. /I 



1-^ ^Jl ^ 



(116) 
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has the Taylor development 



2 °° 1 1 1 1 1 

- + y( - ^ - - - 2Cf ^) . (117) 

5=0 

Again the terms z'^,z,z^ and cancel and shifting each sum by the appro- 



/\gam tne terms z , z, z ana z 
priate amount, eq. (117) writes 

-C'lt'-C'lf-2C'A^ = y^C''^^. (118) 

q=0 ^ 



;+6 
■g+3 



Thus, the third order in perturbations is 

rff = gf-^^C',' + 9192 MCf + 9.^^Cf . (119) 
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